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In this paper, we study the ground state of a one-dimensional exactly solvable model with a 
spiral order. While the model's energy spectra is the same as the one-dimensional transverse field 
Ising model, its ground state manifests spiral order with various periods. The quantum phase 
transition from a spiral-order phase to a paramagnetic phase is investigated in perspectives of 
quantum information science and mechanics. We show that the modes of the ground-state fidelity 
and its susceptibility can tell the change of periodicity around the critical point. We study also 
the spin torsion modulus which defines the coefficient of the potential energy stored under a small 
rotation. We find that at the critical point, it is a constant; while away from the critical point, the 
spin torsion modulus tends to zero. 
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I. INTRODUCTION 

Spiral (or helical) order is a very important phe- 
nomenon in condensed matter physics, especially in soft 
matters. A typical case is the well-known deoxyribonu- 
cleic acid DNA. In solid state, the spiral order can also be 
found in spin systems with competing exchange interac- 
tions. Studies on the spiral order in these systems could 
be traced back to the later of 1950 's[li. Up to now, 
many works on the spiral order have been published 043 ■ 
Among these studies, an interesting example is the one- 
dimensional XXZ spin chain with unparallel boundary 
fields 0- The model's ground state shows a spiral order. 
Moreover, it has been solved exactly using the Bethe- 
ansatz method Q, hence many physical properties can 
be well studied without any approximation. 

However, as far as we know, no one addressed the crit- 
ical phenomena occurring in spin systems with spiral or- 
der in perspective of quantum information science 0. 
So the main motivation of the present paper is to study 
the spiral order transition in terms of the concept of fi- 
delity from the quantum information science. For this 
purpose, we consider a quantum spin chain whose ground 
state consists of a spiral order phase and a paramagnetic 
phase. The quantum phase transition occurring between 
the two phases is studied in perspectives of quantum in- 
formation science and spin torsion modulus. We show 
that the modes of the ground-state fidelity and its sus- 
ceptibility can tell the change of the period of the spiral 
order around the critical point. We study also the spin 
torsion modulus, which defines the coefficient of the po- 
tential energy stored under a small rotation, around the 
critical point. 

The paper is organized as follows: In section HH we 
introduce the model and solve it exactly in a standard 
procedure. The quantum phase transitions from a spiral 
order phase to a paramagnetic phase is also discussed. 



In section IIII1 we study the quantum phase transition 
in term of mode fidelity and its susceptibility. We show 
that the mode fidelity and its susceptibility can not only 
detect the critical point, but also indicate which mode un- 
dergoes a significant change at the critical point. In sec- 
tion [TV] we study the spin torsion modulus of the model 
at the zero temperature. We find away from the critical 
point, the normalized spin torsion modulus tends to zero 
as the system size increases; while at the critical point, it 
is a constant. Finally, a summary is given in section fVl 



II. THE MODEL HAMILTONIAN AND ITS 
EXACT SOLUTION 

We consider a spin chain described by the Hamiltonian 
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where Sj is the twisted 1/2 spin operators. They are 
defined as 
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where <f> is an angle and a = (a x ,a y ,a z ) are the Pauli 
matrices. It is easy to check that s = (s x , s v , s z ) satisfies 
the su(2) Lie algebra, i.e. 
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If we introduce 
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then the Hamiltonian becomes 
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+ (e^+^af4 +1 + h.c)} -h^- (5) 

i=i 

The Hamiltonian either exchanges the state of a pair of 
anti-parallel spins, or flip two upward spins to downward 
or vice versa., so if we define a parity operator 
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and the Hamiltonian cannot change the parity of the 
state, that is. 



\H, P] = 0. 



(7) 



Therefore, we have two subspaces corresponding to parity 
P = ±1 respectively. 

If 4> — 0, the Hamiltonian describes the one- 
dimensional transverse field Ising model with fer- 
romagnetic order 
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While if (j> = 7r, it describes the same model but with 
antiferromagnetic order. Mathematically, if we introduce 
a unitary transformation 



= J] exp 
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The Hamiltonian (TTJ) can be transformed to the one- 
dimensional transverse field Ising model, i.e. 
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Therefore, the energy spectra of the Hamiltonian (TTJ) is 
the same as the quantum Ising model if the corresponding 
boundary conditions are satisfied. 

Nevertheless, what we are interested in this paper is 
the change in the structure of the ground-state wave- 
function around the quantum critical point. For this 
purpose, we still want to take the standard procedure 
to diagonalize the Hamiltonian (JIJ and get the ground- 
state wavefunction explicitly. The procedure consists of 
three transformations, as shown below. 

The Jordan- Wigner transformation: The Jordan- 
Wigner transformation maps 1/2 spins to spinless 
fermions, that is 
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FIG. 1: (color online) A sketch of the spiral order. LEFT: a 
spiral with (f> = 7r/2 and period S = 4. RIGHT: a spiral order 
with cj> = 7T and period S = 2, which actually corresponds to 
an antiferromagnetic order. 
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Then the Hamiltonian becomes 
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If the twist angle satisfies <j> = 2nir/N, n = 0, 1, N — 1, 
the fermionic version of the Hamiltonian has antipcriodic 
boundary conditions for P = 1 and periodic boundary 
conditions for P = —1, respectively. In this case, the 
energy spectra of the Hamiltonian is the same as the 
quantum Ising model. 

The Fourier transformations: Clearly, if h is infinite, 
all spins are aligned to the z direction, then P = 1. So 
the ground state locates in the subspace of P = 1 . Under 
this condition, we can make Fourier transformations 
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e~ ikj c k , 
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where the momentum fcs are chosen as 

(2n+l)7r 
k = . 

N 

The Hamiltonian then becomes 

H = [(2 cos(fc-0) -2h)c\c k 
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-zsin(fc - 0) [c\&_ k+2(t) + c k c- k+ 24,) + h 

The Bogoliubov transformation: The above quadratic 
Hamiltonian can be further diagonalized under the Bo- 
goliubov transformation: 

c k = u k b k + iv k b [ _ k+2c/) , 
4 = u kb\ - iv k b^k+2<t>, 

cl/c+20 = u -k+24,bl k+24> - iv-k+^bk, (16) 

where b k and b\. are also fermionic operator and satisfy 
the same anticommutation relation as c k and c k . Because 
of this, one can find the coefficients in the transformation 
(fl"6|) should satisfy the following condition 

U k = U-k+2tf>, 

v k = u\+vl = l. (17) 

So we can introduce trigonal relation 

v k = sin6» fe , u k = cos9 k . (18) 

Inserting the Bogoliubov transformation into Eq. (fl"5|) . 
the Hamiltonian can be diagonalized under the condi- 
tions 

h - cos(k — 6) 

v / l-2hcos(k-(t)) + h 2 

■ — sin(fc — 6) , . 

sin26» fc = v v ' 19 

v /l-2/icos(fc-0) + /i 2 

Then the Hamiltonian becomes 

H = ^2e(k)(2blb k -l), (20) 

k 

where e(fc) is the dispersion relation 

e(k) = y/l-2hcos(k~(j)) + h 2 . (21) 

From Eq. (|2Tj|) . we can judge that the ground state is a 
vacuum state of b k . Then the ground-state energy can 
be calculated as 

E {^>) = -^2 V 1 - 2/ icos(fc-0) + /i 2 . (22) 

k 

From the dispersion of quasi particles [Eq. (f2~Tj) ]. we 
can see that the system is gapless only at h c = 1 at 
which a quantum phase transition is expected to occur. 
If h 7^ 1, the system is gapped. The two phases can be un- 
derstood from the corresponding limiting cases. If h S> 1, 
all spin are oriented along the z-direction. It is a para- 
magnetic phase. If h ~ 0, the spin chain shows a spiral 




FIG. 2: (color online) 3D surface figure for the mode fidelity 
susceptibility as a function of k and h. Here <f> — 0, N = 40. 

order. The period of the spiral order is determined by <f>. 
In Fig. 1, we show a sketch of spiral orders with cf> = ir/2 
and <f> = 7T respectively. The <j) — n mode actually corre- 
sponds to the antiferromagnetic order. The competition 
between spiral orders and paramagnetic state leads to the 
quantum phase transitions at h c = 1. 

Since the ground state \^o) must satisfy b k \^o) = 0, to 
obtain the ground state, we need to express the operator 
b k {b- k+ 24,) in terms of c k (c k ), 

b k = cos 8 k c k - i sin 9 k c'_ k+2i)l , 
b-k+20 = cos 8 k c- fc+2 + i sin & k c\. (23) 
The condition b k |\l/o) = gives 

a = cos6 k ,b = Q,c = 0,d = is'm8 k . (24) 
The ground state takes the form 
|*oW>= XI (cos9k\0)k\0)- k+ 2^ + ismd k \l) 

A: pairs 

Since the parity here is even, the excited states can 
be obtained by applying 2n b k s with different fcs to the 
ground state. The number of states which belong to this 
subspace is not 2^, but Y^2j ~ 2 JV_1 . With a similar 
procedure, we can obtained another 2 N ~ 1 in the subspace 
of P = —1. Then the total number of state is 2 N and the 
Hilberb space is complete. 

III. MODE FIDELITY, ITS SUSCEPTIBILITY, 
AND SPIRAL ORDER TRANSITIONS 

In this section, we are going to study the quantum 
phase transition occurring at the critical point h c = 1 
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FIG. 3: (color online) 3D surface figure for the mode fidelity 
susceptibility as a function of k and h. Here <j> = tt/2, N = 40. 



FIG. 4: (color online) 3D surface figure for the mode fidelity 
as a function of k and h. Here <f> = tt/2, N = 40, and starting 
h = 0.5. 



in terms of a concept called fidelity emerging from quan- 
tum information science. The fidelity measures the sim- 
ilarity between two states. It has been used to study 
quantum phase transitions in recent years because it can 
describe the structure change of the ground-state wave- 
function as the system varies across a quantum critical 
point fl2l - fl9l |. According to previous publication record, 
it is commonly believed that the fidelity can not only 
detect conventional phase transitions described by Lan- 
dau's symmetry-breaking theorem, but also some uncon- 
ventional phase transitions, such as topological quantum 
phase transitions and Kosterlitz-Thouless phase transi- 
tions. 

For the present model, the fidelity between h and ft/ 
can be calculated as 



F(h,h') = \{* (h')\Mh))\= LI 



cos(6/ fe 
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If 8h = h — ft' is small, F(h,h') shows a minimum around 
the critical point due to the big difference between two 
states coming from the two phases. However, F(h, ft.') 
still depends on 5h, which is rather artificial. In order to 
see the leading contribution to the fidelity, people usually 
expand the fidelity in terms of Sh, i.e. 



F(h,ti) ~ 1 
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XF(h), 
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where xf(^) is the fidelity susceptibility. From Eq. ([26| . 
the fidelity susceptibility can be calculated as 
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Clearly, the fidelity susceptibility takes the same form as 
that of the one-dimensional transverse field Ising model 
except a phase shift in k space. The phase shift does not 
affect the value of the fidelity susceptibility. So around 
the critical point, the fidelity susceptibility has the same 
singular behavior as that in the Ising model [l2j], i-e. 
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If we introduce a, fj,, d, v as the critical exponent of 
XF(h), scaling dimension of XfQi) at the critical point, 
scaling dimension of XF(h) away from the critical point, 
and the critical exponent of the correlation length, they 
satisfy a scaling relation a = (fi — d)v [l6|, 0|. So the 
fidelity susceptibility can not only witness the quantum 
phase transition, but also the universality class of the 
transition. Except for these functions, however, it seems 
that the fidelity susceptibility can not tell more about 
the quantum phase transition. In order to see more de- 
tails occurring around the critical point, we decompose 
the fidelity susceptibility and introduce the mode fidelity 
susceptibility XF(h,k), which defines the changing rate 
of a specific mode as the driving parameter h varies. For 
the present model, the mode fidelity susceptibility takes 
the form 
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XF(h,k) = - 



sin 2 (fc - 



8[l-2/icos(fc-0)+/i 2 ] 5 



(31) 



To see the validity of the mode fidelity susceptibil- 
ity, we firstly let 4> = 0, which corresponds to the one- 
dimensional transverse field Ising model with ferromag- 
netic order. Around the critical point, h c = 1, a quan- 
tum phase transition from a ferromagnetic phase to a 
paramagnetic phase occurs. Since the ferromagnetic or- 
der refers to k = mode, we expect the k = mode 
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FIG. 5: (color online) 3D surface figure for the first derivative 
of mode fidelity as a function of k and h. Here <fr = tt/2, N = 
40, and starting h = 0.5. 
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TABLE I: The critical exponent a, the scaling dimension at 
the critical point /i, and the critical exponent of the correla- 
tion length for the fidelity susceptibility, the relevant mode 
fidelity susceptibility, the fidelity per site, and the mode fi- 
delity. 
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critical point, while other irrelevant 
significant changes. We take a sam- 
as an example, and show the above 
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= N/n, and k = (2m + l)ir/N, then 
l]ir/N . If m — n, we can obtain 



XF(h = h c ,k) , , 



AT2 

8t? 



(32) 



Actually, iff k — <fi ~ Xf(1i>, k) diverges in order of 

N 2 at the critical point. On the other hand, if we let 
k = <ft directly, we find 
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around the critical point. Therefore, for the mode fidelity 
susceptibility, a = 2, /i = 2, d = 0, v = 1, The scal- 
ing relation a — v{[i — d) is satisfied. Clearly XF(h,k) 



diverges more quickly than the normalized fidelity sus- 
ceptibility xfW/N whose critical exponent is 1. On 
the other hand, if fc — is not proportional to 0(1/N), 
XF(h,k) is intensive and independent of the system size 
N. So XF(h, k) shows no singular behavior around the 
critical point. 

As a numerical demonstration, we take a sample of 
N = 40 sites and 4> = tt/2 as an example. We show the 
mode fidelity susceptibility as a function of <f> and h in 
Fig. |31 From the figure, we can see that the mode fi- 
delity susceptibility shows singular behavior at the mode 
4> = 7r/2, and no singular behavior elsewhere. From this 
point of view, the mode fidelity susceptibility might be 
more efficient to describe the quantum phase transition 
occurring in the present model. 

Now we compare the mode fidelity with another re- 
lated and well studied concept, i.e., the fidelity per site 
[l8l |. For the present model, the fidelity per site can be 
calculated as 



\nT(h,ti) = ^lim^i ^ ln[cos(0 fc - (34) 
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It is not difficult to show that the first derivative of the 
fidelity per site [HI 
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around the critical point. Therefore, for the fidelity per 
site, a — 0. 

However, if we focus on the mode of fidelity, 



F{h,h',k) = cos(9 k - e' k ), 



(37) 



it can show more detail about the change in each mode. 
We show F(h,h',k) as a function of hi and k in Fig. @] 
for the case of A^ = 40 and 4> — tt/2. From the figure, we 
can see the k — n/2 decreases the most quickly around 
the critical point. To see the changing rate of the mode 
fidelity, we can take its first derivative 
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Fig. [S] shows a 3D map of dF(h, h! , k)/dh' as a function 
of h' and k. Around the critical point, the most relevant 
mode scales like 
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FIG. 6: (color online) Spin torsion modulus as a function of 
h for various system size. The inset denotes the normalized 
spin system modulus. 



Therefore, for the mode fidelity at k = cf>, a = 0, pt = 
1, v = 1, d = 0, the scaling relation is also satisfied. While 
if k 7^ (f>, dF(h,h',k)/dh is trivial around the critical 
point. That is the mode fidelity can not only tell us which 
mode of the wavefunction undergoes a significant change 
around the critical point, but might be more singular 
than the fidelity per site. 

As a brief summary, we list all the relevant exponents 
in Table U for comparisons. From the table, we can see 
that the mode fidelity and its susceptibility are much 
singular than the fidelity per site and normalized fidelity 
susceptibility. Furthermore, they can tell us which mode 
undergoes a significant change around the critical point. 
Therefore, they are more efficient to describe the critical 
phenomena occurring in this model. 



IV. SPIN TORSION MODULUS 

In this section, we are going to study the mechanical 
response of the spin chain under a small rotation. In 
mechanics, a flexible elastic spring stores mechanical en- 
ergy when it is twisted. As long as the spring is not 
twisted beyond its elastic limit, it obeys an angular form 
of Hooke's law 



-K(j), 



(41) 



where r is the torque exerted by the spring and k is a 
constant called the spring's torsion elastic modulus. The 
potential energy stored in the spring is defined as 



U 



-K0 2 . 



(42) 



For the present model, if we let the spin chain along 
the z-direction, the rotation performed on the spin chain 



behaves like spin torsion. Then we can introduce the spin 
torsion modulus as 



d 2 E 



(43) 
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For arbitray 4>, the energy spectra of the Hamiltonian is 
not the same as the one-dimensional tranverse field Ising 
model. The boundary conditions in Eq. (|T2"j) arc twisted 
boundary conditions. Nevertheless, the Hamiltonian can 
be solved and ground-state energy can be still calculated 
as 



Ea{<f>) 



V 1 - 2hcos(k - (f) + h 2 



(44) 



From Eq. (|43|). the spin torsion modulus of the ground 
state is 
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h 2 [l + cos 2 (fc)] - h(l + h 2 ) cos(fc) 



[1 + h 2 - 2h cos(fc) 



3/2 
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Fig. [5] shows the spin torsion modulus as a function of 
h for various system sizes. From the figure, we can see 
that the normalized modulus is non-zero only at the criti- 
cal point. Away from the critical point, it is almost zero. 
The reason is that in the non-critical region, the sys- 
tems is gapped, the ground state keeps unchanged under 
a small rotation due to the gap protection. Hence only 
at the critical point, the normalized modulus is nonzero. 
This observation can also be obtained exactly. At the 
critical point 



\- 1 



cos 2 (fc) — 2 cos(fc) 



[2 -2cos(fe) 



3/2 



(46) 



which at the infinite limit k/N — 1/tt. While for h ^ 1, 
k/N^0. 

The spin torsion modulus will become the spin stiffness 
if the z-component of the Hamiltonian is conserved. This 
situation is happened at the isotropic case of the one- 
dimensional XY model, whose Hamiltonian reads 



- n 



(47) 



Under the rotation, the XY model can be transformed to 
a fcrmionic model with twisted boundary conditions 



JV 

The spin stiffness can be calculated as(2C 

1 d 2 E{cj)) 



D r . = 



N d(j) 2 



(48) 
(49) 
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For the XY model with twisted boundary conditions, the 
ground state energy is E = — J2k [ cos (k — 0) ~ h]. The 
charge stiffness for h = is 

D c = -. (50) 



V. SUMMARY 

In summary, we studied the quantum phase transition 
occurring in the ground state of a quantum spin chain 
spin with spiral orders in perspectives of quantum infor- 
mation science and mechanical response. We show that 
the mode of fidelity and its susceptibility can not only 
detect the quantum critical point but also indicate which 
modes undergoes a significant change around the transi- 
tion point. Our results manifest that the mode fidelity 
and its susceptibility are much singular that the global 
fidelity and its susceptibility. From this point of view, we 



guess they are more suitable to detect high-order tran- 
sition. We study also the spin torsion modulus of the 
model. The normalized modulus tends to zero in the 
non-critical region, only at the critical point, it becomes 
a constant. Therefore, the spin torsion modulus may 
help us to understand quantum critical phenomena from 
a mechanical point of view. 

Note added. During finishing the work, we received 
a preprint from P. D. Sacramento who used the mode 
fidelity and susceptibility to study a superconducting 
model[2lj]. 
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